Novel Spin Dynamics in a Josephson Junction 
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We address the dynamics of a single spin embedded in the tunneling barrier between two supercon- 
ductors. As a consequence of pair correlations in the superconducting state, the spin displays a rich 
and unusual dynamics. To properly describe the time evolution of the spin we find the generalized 
Wess-Zumino-Witten-Novikov term in the effective action for the spin on the Keldysh contour. The 
superconducting correlations lead to an effective spin action which is non-local in time leading to 
unconventional precessions. Our predictions might be directly tested for macroscopic spin clusters. 
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There is a growing interest in a number of techniques 
that allow one to detect and manipulate a single spin in 
the solid state. A partial list includes optical detection 
of electron spin resonance (ESR) in a single molecule Q , 
tunneling through a quantum dot , and, more recently, 
ESR-scanning tunneling microscopy (ESR-STM) tech- 
nique HQ- Interest in ESR-STM lies in the possibility 
to detect and manipulate a single spin 0,0, which is cru- 
cial in spintronics and quantum information processing. 
In a similar investigation much work has been done to ad- 
dress the coupling, feedback effects, and decoherence in a 
coupled clcctronic-vibrational systems, such as nanomc- 
chanical oscillators and local vibrational modes 0- In a 
previous publication, two of us studied the effect of a pre- 
cessing spin on the supcrcurrent @|. In the present Let- 
ter, we complement this earlier study and focus on how 
a spin is affected by a supercurrent. Our major finding 
is that the spin dynamics is no longer that of simple pre- 
cession. We provide a direct prediction of the expected 
spin dynamics. This unusual spin dynamics character- 
ized by longitudinal oscillations is caused by coupling to 
a Josephson current. Due to spin normalization, these 
novel spin precessions do not lead to any changes in the 
Josephson current. Keldysh contour calculations illus- 
trate that a non-local in time single fermion action is 
also found in situations wherein the single spin cluster is 
replaced by an Anderson impurity 0- As well known, 
in the limit of small hopping amplitudes to and from an 
Anderson impurity, the impurity attains a Kondo like 
character much like that of the single spin which is the 
focus of our attention. Here we consider the origin of 
this rather generic non-locality in time present in the dy- 
namics of a Josephson junction. Our primary focus will 
be on larger spin clusters for which direct measurements 
can be made to probe this novel spin dynamics. 

The model system under consideration is illustrated 
in Fig. ^ R consists of two ideal superconducting leads 
coupled to each other by a single magnetic spin. In the 
presence of a magnetic field, the spin precesses around 
the field direction. We neglect the interaction of the 
spin with two superconducting leads. The Hamiltonian 
for the Josephson junction can be written as [lfil Hll |: 




FIG. 1: Magnetic spin coupled to two superconducting leads. 
In the presence of a magnetic field B, the spin precesses 
around the field direction. As we illustrate here, in a Joseph- 
son junction the standard precession is accompanied by a po- 
lar motion leading to a full blown rigid body like dynamics. 



H = Hl + Hr + Ht ■ The first two terms are the 
Hamiltonians for electrons in the left and right super- 



conducting leads Hl(r) 
3Ew,w[A^(fc(p))cL 
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12k(p);<7 e *(p) c fc(p),<r C fc(p)>' 

2 ^k(p);a,<Tn^™'\*\vjj^ ha cl k( _ p) a ,+K.c.] , where we 
have denoted the electron creation (annihilation) opera- 
tors in the left (L) lead by ct {cka) while those in the 
right (R) lead by (c pr7 ). The quantities k (p) are mo- 
menta, a the spin index, and €k{p),o, A CT(T < (/c(p)) are, re- 
spectively, the single particle energies of conduction elec- 
trons, and the pair potential in the leads. In this work, we 
consider s-wave pairing symmetry in the superconducting 
leads. The two leads are weakly coupled via the tunneling 

Hamiltonian: H T = J2k,p;<T,a'[ T ^'( k 'P) c ka c p^' + H.c] , 
where the matrix element T aa >(k,p) transfer electrons 
through an insulating barrier. When a spin is embedded 
in the tunneling barrier, the tunneling matrix becomes 
a spin operator 0: T = TqS^ + TiS • er CTCT ' , where T 
is a spin-independent tunneling matrix element and 7\ 
is a spin-dependent matrix element originating from the 
direct exchange coupling J of the conduction electron 
to the localized spin S. We take both to be momen- 
tum independent. This is not a crucial assumption and 
is merely introduced to simplify notations. Typically, 
from the expansion of the work function for tunneling, 
ijr ~ J/U, where U is a spin-independent tunneling bar- 
rier || . We further allow a weak external magnetic field 
B ~ 10 2 Gauss. Such a small field will not influence the 
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superconducting state and we may ignore its effect on 
the leads. The Josephson junction with the spin has two 
time scales: (i) The Larmor precession frequency of the 
spin u>l = 9^bB, where g, [Ib are the gyromagnctic ratio 
and Bohr magneton of the conduction electron, respec- 
tively, (ii) The frequency luj = 2eV characterizing the 
Josephson effect when an external voltage V is applied. 

We now derive the effective action via the Keldysh 
technique. If all external fields are the same on both 
forward and backward branches of the Keldysh contour 
(C) then Z = TrT c cxp[-i § c dtH T (t)] = 1, where the 
trace is over both the electron and the spin degrees of 
freedom. We first take a partial trace in Z over the lead 
fermions (the bath) to obtain an effective spin action. 
The Josephson contribution to the resulting spin action 
reads iSS = -\ § c dt § c dt'(TcH T (S{t),t)H T (S(t'),t')), 
much in the spirit of Refs. [12j. 



For brevity, we set A a 



k,p c ka L -P° 



The tunnel- 



ing Hamiltonian of a phase (voltage) biased junction 

Ht = [ToSaa' + ?\S • a aa i\ 

x {A aa , exp(i</>/2) + A\ a , exp(-#/2)) . (1) 

In the presence of a dc voltage bias, <j> = 2eVt. As <j> is 
treated classically (i.e. 4> is the same on the upper and the 
lower branches of the Keldysh contour), the contribution 
tx Tq to SS vanishes. The mixed contribution cx T$T\ 
vanishes due to the singlet spin structure of the s-wave 
superconductor. The only surviving contribution reads 



iSS = -± § c dt § c dt 1 [S(t) • trap] [S(f) • er Sl ] x 
( K {T c A a p{t)A Sl {t'))e^ m ^ 1 + (A,cj>—> A\ -</>)) (2) 

where we keep only the Josephson (off-diagonal) terms. 
The spin structure simplifies for the s-wave case: 



iSS = Tf I dtl dt 1 [S(t) • S(t')] [iD(t,t')\ 
Jc Jc 



(T c A n (t)A l me i:m ^ + (A,, 



(3) 



where iD(t, t') = 

A\ -(/>). Next, we perform the standard Keldysh manip- 
ulations, defining upper and lower spin fields S u,i residing 
on the forward/backward contours and reducing the time 
ordered integral over Keldysh contour to the integral over 
forward running time at the cost of making the Green's 
function G a 2 x 2 matrix. Finally, after a rotation to the 
classical and quantum components 

Si = (S'" + S')/2, S 2 = S u -S l , Si-S 2 = 0, (4) 

we obtain iSS = iSn + Si, where 

iS R = J J dtdt'S 2 (t)-S 1 {t')[iK 12 (t,t')}. (5) 
with K 12 {t,t') = If [D R (t,t') + D A (t',t)}, and 

Sl = \j J dtdt ' S2{ $ ■ Mt%K22(t,t% (6) 




FIG. 2: The unit sphere for the vectors n°''(i) is shown. The 
contour C is the Keldysh contour for the forward (u) and 
backward (I) evolution. To properly describe the spin on this 
closed contour we analyze the WZ WN term, see Eq. 113H . For 
clarity, we draw a small piece of the closed trajectories. 



with K 22 (t,t') = T?D K {t,t'). The retarded (R) and 
Keldysh (K) components are defined via 



iD R (t,t') = Q(t-t') 



x([A n (t), AuCf)]-)***^ -cc. (7) 



iD K (t,t') = ({A(t) n ,A u (t')} + y 



+ c.c. (8) 



The advanced component in our case is simply 
iD A (t,t') = iD R {t',t). The kernels K 12 and K 22 are 
readily calculated at T = from 



(A(t) n A u (t'))=J2 



k . p 



4EkE p 



-i(E k +E„)(t-t') 



(9) 



where E k{p) = y/(e k ( p ) - E F ) 2 + |A| 2 . As only frequen- 
cies higher than 2 A are present in Si (Eq. (JSJ) van- 
ishes for slow fluctuations of S2 unlike the real part of the 
action Sr (due to the presence of 8 function in Eq. 

To properly describe the dynamics of the spin fields Si 
and S2, we em ploy the path integral representation for 
the spin fields [13|. The action for a free spin consists 
of two terms Sq = g^B § c dtYi ■ S + Swzwn- The sec- 
ond, Wcss-Zumino-Wittcn-Novikov (WZWN), term de- 
scribes the Berry phase accumulated by the spin as a 
result of motion of the spin on the sphere We gen- 
eralize this action for nonequlibrium dynamics within 
the Keldysh contour formalism (Fig. |2Jl. We write the 
WZWN term on both forward and backward contours, 
and consequently analyze matters in terms of Si and S2, 

iSwzwN = / dT J dt[S u (t,r) ■ (d T S u (t,r) 

xd t S u (t,T))-(u^l)] . (10) 

The relative minus sign stems from the backward time 
ordering on the return part of C. The spins S u ' 1 = Sn u ' 1 
with S the magnitude of the spin and n u ' 1 a unit vec- 
tor field. The additional integral over r permits us to 
express the action in a local form. At r = we set 
the spin to point along the z direction at all times; at 
r = I the spin field corresponds to the physical configu- 
rations. Each of the individual WZWN phases (for both 
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the forward (u) and backward (I) branches) is the spin 
magnitude (S) multiplied by the areas spanned by the 
trajectories n u,l (t) on the unit sphere. The WZWN term 
contains odd powers of S2. Insofar as the WZWN term 
of Ea. H10|) is concerned, the standard Keldysh transfor- 
mation to the two classical and quantum fields, Si and 
S2, mirrors the decomposition of the spin in an antifer- 
romagnet (AF) to the two orthogonal slow and fast fields 
[lij . The difference between the two individual WZWN 
terms in Eq. (|10J) is the area spanned between the forward 
and backward time trajectories. The magnitude of this 
area traced between times t and t + St is 



WZWN 



= S(5t)5n • (n x <9 t n). 



(11) 



Here, the variation between the forward and backward 
trajectory at a given instant of time is <5n = S2/S. In 
Ea.l|ll[l. wc note that for small variations between the 
forward and backward trajectories, n = Si/S. The 
WZWN action on the Keldysh loop may be expressed 
as 



Swzwn — ^2 / 



(Si x 9 f Si). 



(12) 



The real part of the total action, S c i = Sq + Sr, de- 
termines the quasi-classical equation of motion. The 
imaginary part of the action Si is usually responsible for 
the Langevin stochastic term. In our case, however, the 
Langevin term is suppressed at frequencies much lower 
than A. The fluctuations of S2 are not suppressed by 
the bath and no dissipation appears. The only reason to 
consider small fluctuations of S2 is the dominance of the 
quasi-classical trajectories determined by S c i. Thus our 
analysis applies for large spins. The classical action 

Sd = Swzwn + g^B J dtH ■ S2 
+ j dt J dt'K 12 {t,t')S 2 (t) • Si(i'). (13) 

As the spin dynamics is much slower as compared to elec- 
tronic processes, we set Si(i') ~ Si(t) + (t' — t)dSi/dt. 
The variational equations 6S c i/SS2(t) = imply 



dn dn . 

an x — — sin a; jt 



dt 



dt 



(14) 



and a -- 
-eV)- 2 } 



where henceforth we denote Si by S = S11, 

s E k , P + e p ~ eV ) ~ 2 ~ ( E * + e p- 

The non-dissipative term proportional to a in Eq. (|14|> 
arises from superconducting retardations. Both Sx| 
and I ~ sin ujjt are odd in time and their product is al- 
lowed in the equation of motion. Due to the spectral gap, 
dissipation is faint at low temperatures and frequencies. 

The classical equation of motion, Eq. (|14|) . implies 
that n ■ dn/dt = 0, as it must to be consistent with 
the parameterization of the spin on the sphere, Sri = 
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FIG. 3: The spin dynamics. Top panels: the azimuthal ((f>) 
and polar (8) angles as a function of time for an exaggerated 
a — 0.4 (red) versus the standard precessions in the absence 
of the coupling (a = in blue). Note that unlike the stan- 
dard case, the semi-classical spin trajectory is not confined 
to planar motions. The polar angle oscillates between its two 
extreme values 81 and 82. In the bottom panels we display 
(S x (t)) and (S z (t)). Here we take ujl = uj. 



S(s'm 6 cos (j>, sin 9 sin 4>, cos 0) . Orienting the z-axis along 
the external magnetic held B, Eq. (|14(l reads 



d(f> 
~dt 
dO 

~dt 



1 + a 2 sin 2 (iojt) ' 

^ ■ a ■ * 
-a — sin a sin uj it . 

dt 



(15) 



For a spin initially oriented at an angle #0 relative to B, 

m = - 



:tan + a 2 tan(wji)], 



6(t) = 2 tan" 1 



(1 -ccos(wjt))(l + c) 17 



(l + ccos(wjt))(l-c) 



tan 



with c = \a\/ yT+c? and 7 = —au>L/2uijc. For a <C 1, 
4> ~ — u>iX and 6 ~ 6q — a(uiL/uJj) sin#o cosuijt. Canon- 
ically, whenever a single spin is subjected to a uniform 
magnetic held, the spin processes azimuthally with a fre- 
quency u>l. In & Joscphson junction, however, the spin 
exhibits additional polar (9) displacements. The result- 
ing dynamics may be likened to that of a rotating rigid 
top. The Joscphson current leads to a full non-planar gy- 
roscopic motion (nutation) of the spin much like that gen- 
erated by applied "torques" on a mechanical top. When 
B = 0, the solutions do not allow for any spin dynam- 
ics. The possibility remains, nonetheless, that coupling 
to random magnetic field (bath) will have a nontrivial 
effect on the spin in presence of a Josephson current. In 
Fig. |3 we display the resulting dynamics for the spin. 
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FIG. 4: The resulting spin motion on the unit sphere in the 
general case. As in the motion of classical spinning top, the 
spin exhibits undulations along the polar direction. 



tant differences with the case of Josephson effect in the 
presence of a vibrational mode [l5| . The normalization 
of the spin disallows resonance effects at lul = ujj. The 
polar oscillations do not lead to additional harmonics in 
the Josephson current. Nevertheless, the dynamics may 
be monitored by SQUID measurements. 

We thank Yu. Makhlin for useful discussions. This 
work was supported by the US DOE. 



A schematic of the generic spin motion is displayed in 
Fig. 0] Similar to a classical spinning top, the spin wob- 
bles along the polar direction in addition to azimuthal 
rotations. Similar dynamics is even expected for a quan- 
tum S = 1/2. In a spin coherent state path integral, the 
Schrodinger equation is essentially classical. 

The Josephson current in the presence of a process- 
ing spin has been discussed in ||. Within our approx- 
imation, the terms in the action that may affect the 
current are proportional to S 2 = 1 and S • ^ = 0. 
Consequently, the Josephson current remains unmodu- 
lated. Nevertheless, we may directly monitor the spin 
dynamics by SQUID measurements. The spin motion 
generates a time dependent magnetic field, <5B(r,t) = 
^■[3r(r • m(t)) — r 2 m(i)]/r 5 , superimposed against the 
constant external field B. Here, r is the position from the 
spin and the magnetic moment m(t) = g[iBSn(t). For a 
ferromagnetic cluster of spin S = 100, a detectable field 
SB ~ 10 -10 Tesla appears in a SQUID placed a micron 
away. For a SQUID loop of micron dimensions, the cor- 
responding flux variation (5$ ~ 10 — 7 <E>q (with $o a flux 
quantum) — within reach of modern SQUIDs. For such 
a setup with (Ti/Tq) ~ lO^ 1 , typical critical Josephson 
current 4° } - 10 fiA, |A| = 1.0 meV, and eV ~ 10" 3 | A|, 
we find that a ~ 0.1. As n x = sin 9 cos <f>, with a sim- 
ilar relation for n y , the spin components orthogonal to 
B vary, to first order in a, with a Fourier component at 
frequency \u>l ± leading to a discernable signal in 
the magnetic field (B + SB). For a field B ~ 200 Gauss, 
Wi ~ 560 MHz wherein for a substantial voltage (wj) 
range, a new side band will appear at \ul — uj\ whose 
magnitude may be tuned to 0(10 — 100) MHz. This mea- 
surable frequency is markedly different from that associ- 
ated with standard Larmor frequency (ujl) precessions. 

Conclusion. By analyzing WZWN phases within the 
non-equilibrium Keldysh framework, we find novel non- 
planar spin dynamics in the presence of a tunneling 
Josephson current. The coupling of the spin to a super- 
conducting bath produces non-damping spin retardation. 
This retardation results from additional non-dissipativc 
terms in the spin equations of motion. There are impor- 



[1] J. Koehler et al, Nature 363, 342 (1993); J. Wrachtrup 
et al, ibid. 363, 244 (1993); Phys. Rev. Lett. 71, 3565 
(1993). 

[2] H.-A. Engel and D. Loss, Phys. Rev. Lett. 86, 4648 

(2001) ; Phys. Rev. B 65, 195321 (2002). 

[3] Y. Manassen et al, Phys. Rev. Lett. 62, 2531 (1989); D. 

Shachal and Y. Manassen, Phys. Rev. B 46, 4795 (1992); 

Y. Manassen, J. Magnetic Reson. 126, 133 (1997); Y. 

Manassen, I. Mukhopadhyay, and N. Ramesh Rao, Phys. 

Rev. B 61, 16223 (2000). 
[4] C. Durkan and M. E. Welland, Appl. Phys. Lett. 80, 458 

(2002) . 

[5] H. Manoharan, Nature 416, 24 (2002); H. Manoharan, 
C. P. Lutz, and D. Eigler, Nature 403, 512 (2000). 

[6] A. V. Balatsky and I. Martin, Quan. Inform. Process. 1, 
53 (2002); A. V. Balatsky, Y. Manassen, and R. Salem, 
Phys. Rev B 66, 195416 (2002). 

[7] D. Mozyrsky and I. Martin, Phys. Rev. Lett. 89, 018301 

(2002) . 

[8] J.-X. Zhu and A. V. Balatsky, Phys. Rev. B 67, 174505 

(2003) . 

[9] Y. Avishai, A. Golub, and A. D. Zaikin, Phys. Rev. B 

63, 134515 (2001). 
[10] V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 

486 (1963); 11, 104(E) (1963). 
[11] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991); 

V. Geshkenbein and A. I. Larkin, JETP Lett. 43, 395 

(1986); D. A. Wollman et. al, Phys. Rev. Lett. 71, 2134 

(1993). 

[12] U. Eckern, G. Schon, and V. Ambegaokar, Phys. Rev. 
B 30, 6419 (1984); A.I. Larkin and Yu. N. Ovchinnikov, 
Phys. Rev. B 28, 6281 (1983). 

[13] E. Fradkin, Field Theories of Condensed Matter Sys- 
tems (Addison- Wesley, Redwood City, 1991); S. Sachdev, 
Quantum Phase Transitions (Cambridge University 
Press, London, 1999). 

[14] These two orthogonal AF fields represent (i) the slowly 
varying staggered spin field (the antiferromagnetic stag- 
gered moment m taking on the role of Si) and (ii) the 
rapidly oscillating uniform spin field 1 (paralleling our 
S2). In the antiferromagnetic correspondence, the two 
forward time spin trajectories at two nearest neighbor 
AF sites become the two forward (u) and backward (I) 
single spin trajectories of the non-equilibrium problem. 
This staggered doubling correspondence is general. 

[15] J.-X. Zhu, Z. Nussinov, and A. V. Balatsky, 
|cond-mat/0306107| 



4 



